Spin-1 chain with spin-1/2 excitations in the bulk 
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We present a spin-1 chain with a Hamiltonian which has three exactly solvable ground states. 
Two of these are fully dimerized, analogous to the Majumdar- Ghosh (MG) states of a spin-1/2 
chain, while the third is of the AfHeck-Kennedy-Lieb-Tasaki (AKLT) type. We use variational and 
numerical methods to study the low-energy excitations which interpolate between these ground 
states in different ways. In particular, there is a spin-1/2 excitation which interpolates between the 
MG and AKLT ground states; this is the lowest excitation of the system and it has a surprisingly 
small gap. We discuss generalizations of our model of spin fractionalization to higher spin chains 
and higher dimensions. 

PACS numbers: 75.10.Jm, 75.10.Pq 
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I. INTRODUCTION 

Quantum spin systems in one dimension have been 
studied extensively for many years. In some seminal pa- 
pers, Haldane predicted theoretically that integer spin 
chains with nearest neighbor Heisenberg antiferromag- 
netic interactions should have a gap between the ground 
state and the first excited stated; this was then ob- 
served experimentally in a spin-1 system^ii^i^ and con- 
firmed numerically^i^i. Haldane's analysis used a field 
theoretic description of the long-distance and low-energy 
modes of the spin system^'^'^"'^^. Affleck, Kennedy, Lieb 
and Tasaki (AKLT) then showed that the ground state 
of the spin-1 chain can be variationally understood as a 
state in which each spin-1 is thought of as a symmetric 
combination of two spin-l/2's, and the two spin-l/2's 
at each site form a singlet with the spin-l/2's of the 
neighboring sites"'^^. The excitations are given by vari- 
ational states in which one of these singlets is replaced 
by a triplet. It was shown later that the AKLT state can 
be written as a matrix product statei^. 

If a spin chain has sufficiently strong next-nearest- 
neighbor interactions, the system is frustrated and its 
low-energy properties can be quite different from those 
of the unfrustrated system. For instance, the spin-1/2 
chain with both nearest neighbor (Ji) and next-nearest 
neighbor (J2) antiferromagnetic interactions is gapless if 
J2 = 0, but is gapped if J2/J1 > 0.241liiii^. In the latter 
case, the ground state is doubly degenerate as expected 
by the Lieb-Schultz-Mattis theorerrki^. In particular, for 
the Majumdar- Ghosh (MG) model given by J2/J1 — 1/2, 
the ground states are exactly solvable^^ and consist of 
products of nearest-neighbor singlet states as will be de- 
scribed below. The lowest excited states then consist of 
spin-l/2's interpolating between the two ground states^^. 
Hence the excitations of the MG model have spin 1/2 in 
contrast to the excitations of the AKLT model which 
have spin 1. 

The excitations described above exist in the bulk; they 
contribute to thermodynamic quantities like the mag- 
netic susceptibility and the specific heat. In addition to 
these excitations, a gapped chain with a finite number of 



sites may also have degrees of freedom localized at the 
edges. For instance, the AKLT model on an open chain 
has spin-1/2 degrees of freedom at the edges^; these can 
be thought of as remnants of the two spin-l/2's of which 
each spin-1 is composed. These edge degrees of freedom 
have been studied using field theoretic methods^^. It 
may be interesting to consider spin-1 chains which have 
spin-1/2 excitations in the bulk, as this would provide an 
example of spin fractionalization. 

Spin fractionalization was first proposed by Fad- 
deev and Takhtajan in the spin-1/2 antiferromagnetic 
chain2£; the idea is that the elementary excitations, 
called spinous, carry spin-1/2. This was confirmed 
experimentally in a one-dimensional spin-1/2 system 
KCuF^^^. It was later shown by Anderson and others 
that spin fractionalization can also occur in higher di- 
mensional systems with resonating valence bond ground 
states2Si2^. This idea has been used to understand the 
low- lying excitations in a two-dimensional spin-1/2 sys- 
tem C S2CuCljr^^. In contrast to these examples of spin 
fractionalization in spin-1/2 systems, we are proposing a 
model of spin fractionalization in higher spin systems in 
this paper. 

A spin-1/2 excitation existing in the bulk of a spin-1 
chain must clearly have two different ground states on its 
left and right. For instance, the ground state on the left 
could be of the MG type in which each spin-1 forms a 
singlet with one of its neighbors, while the ground state 
on the right could be of the AKLT type. The spin-1/2 
excitation can then be thought of as the edge degree of 
freedom of the AKLT part of the chain. To realize this 
kind of an excitation, we require a Hamiltonian for which 
both MG and AKLT states are ground states. We will 
present such a Hamiltonian in Sec. II; it contains inter- 
actions involving three neighboring sites. We will present 
a variational estimate of different possible excitations of 
the model, and will show that the spin-1/2 excitation has 
the lowest variational energy. In Sec. Ill, we will present 
numerical results for finite chains, with both open and 
periodic boundary conditions. These will confirm that 
the spin-1/2 excitations indeed have the lowest energy; 
with periodic boundary conditions, such excitations must 
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occur in pairs. In Sec. IV, we will discuss how our model 
can be generalized to higher spins and higher dimensions, 
i.e., how one can construct models which have spin S at 
each site and spin S' excitations in the bulk, with S' < S. 
We will make some concluding remarks in Sec. V. 



and |//) = Yi \Si2n-l,2n)) . (2) 



These are generalizations of the two ground states of the 
spin- 1/2 chain at the MG point^. 



II. A SPIN-1 CHAIN 



A. Hamiltonian and ground states 



We will first present what appears to be the simplest 
Hamiltonian of an infinite spin-1 chain which has ex- 
actly three ground states. This Hamiltonian is motivated 
by the following arguments. Given three spin-l's 81,82 
and 5*3, let us define the projection operators Ps which 
projects on to states with total spin S, where S can be 
0, 1, 2 or 3. Now consider a three-spin Hamiltonian of 
the form h = C2P2 + C3P3, where 02,03 > 0. The ground 
states of h are all the states whose total spin is equal to 
or 1; all such states have zero energy. All the excited 
states have strictly positive energies. If we think of each 
of the spin-l's as being a triplet combination of two spin- 
1/2 's, these ground states correspond to states in which 
at least four of the six spin-l/2's form singlets amongst 
each other. The remaining two spin-l/2's can at most 
form a total spin of 1, no matter how they combine with 
each other. Now, a particular Hamiltonian of the above 
type is = Sl^{Sl, - 2), where Si, = {S, + S2 + S^)'; 
this corresponds to the coefficients C2 = 24 and C3 — 120. 
This is the simplest Hamiltonian with ground state spins 
being equal to and 1 in the sense that it has the lowest 
possible powers of the spin operators Si. 

We now consider a Hamiltonian for the spin-1 chain of 
the form 

H ^ J 

n 

where /i„ = [Sn-i + Sn + Sn+i)"^ 

X [{Sn-l+S,, + Sn+lf - 2] . (1) 

(We will set the exchange constant J equal to 1). The 
ground states of this Hamiltonian must have at least two 
singlet bonds within every group of three neighboring 
spins. It is then easy to see that there are three degen- 
erate ground states with zero energy of the forms shown 
in Fig. 1. The analytical expressions for these three 
states are as follows. Let us define the singlet combi- 
nation of two spin-l's at sites m and n as \S{m,n)) = 
[|l,-l)mn - |0,0)m„ + | - l,l)m„]/V3, whcrc wc have 
used the 5^ components to label the states. Then the 
first two ground states of ([T]) are given by tensor prod- 
ucts of singlets between nearest neighbors of the form 



|/) = n l'S(2n,2n+l)) 
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FIG. 1: The three degenerate ground states. Each solid circle 
represents a spin 1/2, and the lines denote singlet bonds. 

The third ground state of Eq. (P) is the AKLT state. 
This can be written as a matrix product state^'^. At a 
site n, let us define the matrix 



Then the AKLT state is given by the matrix product 



im) = 



n Mr,. 



(4) 



The matrix in Eq. ^ is motivated as foUowa^^. For 
a spin-1/2 object, we can use u — cos(6'/2)e"'^/^ and 
V = sin(6'/2)e^*"*/^ to describe the spin-up and spin- 
down states respectively. The spin operators are given 
by S^ = {udu — vdv)/2, S*"*" — udv, and S^ = 
the total spin is S* = {udu + vdy)/2. The inner product 
in the (u, v) space is defined by the integration measure 
dn, — sin 9 d9d(j)/ {An). The correctly normalized spin-1/2 
states are given by |l/2) = \/2u and | — 1/2) = \/2v. 
For a spin-1 object, the normalized states are given by 
|1) = V3w^ |0) = VQuv, and | - 1) = VSv^. A singlet 
formed by spin-l/2's at sites n and ri + 1 is given by 



UiV^+l - ViUi+l = (uiV.A [ ) 



(5) 



The matrix in Eq. Q is obtained by combining a column 
and a row for site n as 



M„ V2 I _ I ( Un Vr, 



(6) 



The normalization of M„ has been chosen so that the 
norm of the AKLT state in Eq. ^ is given by 



Tr ( 1/3 =1 in the limit N ^ 00 . (7) 
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The three states defined in Eqs. ([2]) and (|4]) are or- 
thonormal for the infinite chain. We do not have an ana- 
lytical proof that these are the only ground states of Eq. 
([T|). However, we will provide numerical evidence in Sec. 
Ill that there are no other ground states, except for some 
additional degeneracies in open chains due to degrees of 
freedom at the edges. 

The structure factor in a ground state is given by 



S{q) 



(8) 



where N is the number of sites in the chain, and we 
eventually have to take the limit N oo. In the three 
ground states given above, we find that-^^ 



and S^^\q) 



6 (1 - 



= 2(1 
cosq) 



cosq) , 



3 cos (7 



(9) 



B. Excited states 



We will now study the excited states using a varia- 
tional techniquei^i^liS^. Given two ground states A and 
B, which could be any of the states /, // or ///, one 
can consider a 'domain wall' state |n) which interpo- 
lates between the two at site n. We can then superpose 
such states to form momentum eigenstates as shown 
below, and obtain a variational estimate of the energy 
Evar{k) = {k\H\k) / {k\k) . We will now do this for vari- 
ous possible combinations of the two ground states A and 
B on the left and right. There are four different cases to 
consider. In each case, we will form an excited state by 
breaking as few singlet bonds as possible. 



(1) 



(ii) 



(iii) 



(iv) 



FIG. 2: Various possible excitations interpolating between 
different ground states. The lines denote singlet bonds, and 
each isolated circle denotes a free spin 1/2. 

(i) We first consider a state interpolating between ground 
states / on the left and // on the right as shown in Fig. 
2 (i). This is given by 

n 

\2n{I,II)) = II |5(2m-2,2m-l))(^|l)2„ 



Yl |S'(2m + 1,2to + 2)) 



(10) 



This is a state with SL^ = 1. We then find that 



(2m(/,//)|2n(/,//)> = (l/3)l"-"l , 
{2m{I,lT)\H\2n{I,H)) = 40 5™,„ . (11) 



If we form the momentum eigenstate 
|fc) = ^ e**^^" \2n] 



(12) 



we find that 



{k\k) 



2N 



5 - 3cos(2fc) ' 
and {k\H\k) = 20N . (13) 

From Eq. (|13p . the variational energy is given by 

E^ar{k) = 10 [5 - 3cos(2fc)] . (14) 

The minimum of this lies at fc = 0, where Eyar{0) — 20. 
(ii) Next we consider a state interpolating between 
ground states / on both the left and the right as shown 
in Fig. 2 (ii). This is obtained by replacing a singlet 
\S(2n, 2n + 1)) by a triplet. We thus have 

n 

|2n(/,/)) Yl |5(2to- 2,2m- 1)) 



m— — oo 



—7= 0)2n,2n+l ~ |0, l)2n.2n+l] 

v2 



Y[ |S'(2m,2m+l)) 



(15) 



m— n+1 

This is a state with S^^^ = 1. We find that 



{2m{I,I)\2n{I,I)) = 6 

on 

{2m{I,I)\H\2n{I,I)) = — 5.m^n ■ (16) 

A momentum eigenstate defined as in Eq. (fT2l) satisfies 

N AON 
{k\k) = - , and {k\H\k) = ^ . (17) 



Hence the variational energy is 

Evar{k) = 80/3 ~ 26.67 



(18) 



independent of the value of fc. 

(iii) We now consider a state interpolating between 
ground states /// on the left and / on the right as shown 
in Fig. 2 (iii). The ground state /// must end with one 
singlet bond between the spin-l/2's at site 2n and 2n-|-l, 
along with a free spin 1/2 at site 2n + 1. We therefore 
take a state which is of the AKLT type from —oo to site 
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2n; this is followed by a column multiplied by a free spin 
1/2 at the site 2n + 1 of the form 



This is a state with Sf, 



1. We then find that 



The choice of U2n+i, rather than V2n+i, as the free spin 
1/2 at the end of the AKLT region makes this a state 
with S^gf — 1/2. The free spin is then followed on the 
right by the ground state /. The complete state is thus 
given by 



|2n(///,/)) 



n 




Y[ |S'(2m,2m + l)) 



(20) 



m— n+l 



We then find that 



{2m{III,I)\2n{III,I)) = (-l/VS)!""™! 

QQ 

{2m{III,I)\H\2n{III,T}) = — 5™,„ . 

y 



(21) 



A momentum eigenstate defined as in Eq. (fT^ satisfies 

N 



and 



(fc|fc) 
\H\k) 



2 [2 + y3cos(2/c)] ' 
40iV 
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Hence the variational energy is 
80 



[2 + \/3cos(2fc)] 



(22) 



(23) 



The minimum of this lies at A: = 7r/2, where E^ari^^ l"^) — 
2.38. 

(iv) Finally, we consider a state interpolating between 
ground states /// lying on both the left and the right as 
shown in Fig. 2 (iv). We take the AKLT state on the 
left to be of the same form as the one discussed around 
Eq. (fT9|) . with 2n replaced by n — 1. The state on the 
right begins with a free spin 1/2 multiplying a row at site 
n + 1 of the form 



\/2 Un+l ( Un+1 Vn+1 ) 

= ( v/273 (l|„+i 7173 (0|„+1 



(24) 



This is then followed by a state of the AKLT type from 
site n + 2 to 00. The complete state is thus given by 



\n{IIIJII)) 



n 




{m{III,III)\n{III,III)) 

{m{III, III)\H\n{III, III)) 

_ 320 _ ^ ( X 

27 ™'" 27 

A momentum eigenstate defined as 



■^m,n4-l 



satisfies 



\k) 
m 

and {k\H\k) 



E 



(26) 



(27) 



1 
2 

160iV 
27 



- cos k I N 



(2 — cos/c) 



Hence the variational energy is 



Eyar (A) 



320 (2 - cos/c) 
9 (3 - 2cosfc) 



(28) 



(29) 



TT, where Eyar{T^) 



The minimum of this lies at k 
21.33. 

A comparison between the four kinds of excitations dis- 
cussed above shows that the gaps of excitations of type 
(i), (ii) and (iv) are given by 20, 26.67 and 21.33 respec- 
tively, while excitation (iii) has a gap of only 2.38. We 
note that excitation (i) leaves one triangle unsaturated 
by two bonds, i.e., one group of three neighboring spins 
has no singlet bonds within themselves; this can be seen 
in Fig. 2. Excitations (ii) and (iv) both leave two tri- 
angles unsaturated by one bond each. Excitation (iii) 
leaves one triangle unsaturated by one bond. The min- 
imum energy excitation is of type (iii) which represents 
a 'domain wall' interpolating between ground state / (or 
//) and ///, i.e., between ground states of the MG and 
AKLT types. The gap of 2.38 for this state is much less 
than the excitation energy of 24 of the three-state Hamil- 
tonian hn appearing in Eq. ([1]). Further, this state has 
spin 1/2 arising from the free spin 1/2 described around 
Eq. dlSl). 

We have not tried to improve our variational calcu- 
lations by considering more extended states which in- 
terpolate between the different ground states. Such ex- 
tended states do not seem to greatly improve the energy 
estimate^; this is because our ground states have fairly 
short correlation lengths. Further, we will see in Sec. Ill 
that the numerical result for the lowest excitation gap is 
not very different from the variational estimate obtained 
above in Eq. (^5)) . 



III. NUMERICAL RESULTS 



/273 (1 



n+l 



I n+l 



n • (25) 



m— n+2 



We will now study the model defined in Eq. ([T]) us- 
ing exact diagonalization of finite chains, with both open 
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and periodic boundary conditions (PBC). We will check 
whether the three states discussed in Sec. II. A are the 
only ground states, and also what the lowest excitation 
energy is. If the spin-1/2 excitations described in Sec. 
II. B are indeed the lowest energy excitations with a gap 
AE, we would expect the gap for open chains to be given 
by AE while the gap for a chain with PBC should be 
2AE. This is because an open chain may have a single 
spin-1/2 excitation with a gap in the bulk, and a gap- 
less spin-1/2 degree of freedom localized near one of the 
edges which compensates for the spin 1/2 in the bulk. 
But a chain with PBC can only have excitations in the 
bulk which have integer values of S^^^ ; hence these must 
consist of at least two spin-1/2 excitations. 

We have studied chains with N ranging from 5 to 10. 
In the exact diagonalization procedure, we used the quan- 
tum number S^^^ and symmetry under parity to reduce 
the sizes of the Hilbert spaces. For open chains with an 
even number of sites, the degeneracy of ground states is 
found to be 14. This confirms that the three states dis- 
cussed in Sec II. A exhausts the list of all ground states 
since it can be understood as follows using Fig. 1. There 
is one state of type I, 9 states of type II (there are two 
unpaired spin-l's at the edges giving a degeneracy of 3^), 
and 4 states of type III (the two dangling spin-l/2's at the 
edges give a degeneracy of 2^). For an open chain with 
an odd number of sites, we find 10 degenerate ground 
states. This can be counted as 3 states each of types 
I and II arising from an unpaired spin-1 at one of the 
edges, and 4 states of type III due to the two dangling 
spin-l/2's at the edges. 

For chains with PBC and an even number of sites, we 
expect 3 degenerate ground states corresponding to each 
of the three types. For an odd number of sites, ground 
states of types I and II are not allowed because they 
would leave one triangle unsaturated; thus we expect a 
unique ground state of type III. These expectations have 
been confirmed by the numerics. 

Next we consider the first excited state. Fig. 3 shows 
the energy gaps as a function of the chain length N, for 
chains with PBC (upper two lines) and for open chains 
(lower two lines) . Although the results differ significantly 
between even and odd values of N, they extrapolate to 
about the same values for ^ oo. We have fitted the 
gaps to the form AE{N) = AE{oo) + a/N'^. The reason 
for this fitting form is that an excited state with a gap is 
expected to behave like a particle in a box^>^; in a system 
of length N, the leading A^-dependent term in the energy 
of such an object is The inset of Fig. 3 shows the 

ratio Aper / Aopen as a function of for even and odd 
values of N; the lines in the inset are obtained by taking 
the ratio of the fitted lines in the main figure. 

In Table 1, we summarize the results shown in Fig. 
3 by listing the gap for various values of A^ for open 
(Aopen) and periodic (Aper) boundary conditions as well 
as the ratio Aper/ Aopen- We see that the gap for the 
open chain extrapolates to a value of about AE = 2.5 
which is not very different from the value of 2.38 obtained 




0.005 0.01 0.015 0.02 „0.025 0.03 0.035 0.04 0.045 

FIG. 3: Energy gap as a function of where N is the 

chain length. The upper two lines are for periodic bound- 
ary conditions, while the lower two lines are for open chains. 
Dashed and solid lines denote even and odd values of A'^ re- 
spectively. The inset shows the ratio Aper/Aopen as a function 
of l/A^.' 



N 


^open 


^per 


^per / ^open 


5 


6.08 


9.28 


1.53 


6 


6.29 


11.28 


1.79 


7 


4.34 


7.09 


1.63 


8 


4.63 


8.72 


1.88 


9 


3.61 


6.30 


1.75 


10 


3.82 


7.29 


1.91 


^odd 


2.51 


4.93 


1.96 




2.45 


5.14 


2.10 



TABLE I: Gaps for chains with open and periodic boundary 
conditions for different chain lengths. The last two lines give 
the gaps extrapolated to the thermodynamic limit for an odd 
and even number of sites respectively. 



variationally in Eq. (|23p. Further, the gap for the chain 
with PBC extrapolates to a value which is about twice 
that of the open chain gap. This implies, for instance, 
that there is no bound state of two spin-1/2 excitations 
which has an energy which is significantly less than 2AE. 

For open chains, we find that the total spin of the 
lowest excitation is Stot — 1 for even A^ and Stot ~ 2 for 
odd A^. The latter value can be understood as follows: 
If this excitation is the state (iii) discussed in Sec. II. B 
(see Fig. 2 (iii)), which interpolates between AKLT and 
a fully dimerized ground state, then it is possible to have 
an unpaired spin-1 at the edge of the fully dimerized side 
without costing any energy. This edge spin can combine 
with the spin-1/2 at the edge of the AKLT side and the 
spin-1/2 in the bulk to form Stot = 2. 
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IV. GENERALIZATIONS 

We can construct models involving higher spins or 
higher dimensions in which excitations in the bulk can 
carry spins which are a fraction of the spin at each site. 
We will discuss some examples below. 



states with any two values of pi and p2 lying in the 
range < P2 < Pi < 25*. We thus see that the excited 
states of this spin-S' chain can have any value of the spin 
{pi — P2)/2 going from 1/2 to S. 

B. Higher dimensional models 



A. Higher spin chains 

The idea of a Hamiltonian with multiple ground states 
in which there are varying numbers of singlet bonds be- 
tween neighboring sites can be generalized to higher spin 
chains. Consider a chain of spin-S* sites with a Hamilto- 
nian such that all ground states must have at least 25* 
singlet bonds amongst every group of three neighboring 
sites. In analogy with Eq. ([T]), we can write such a 
Hamiltonian as H = hn, where hn is a sum of projec- 
tion operators on to values of total spins ranging from 
S" -f 1 to SiS* for sites n—l,n and n+1. A state in which 
there are p singlet bonds between sites 2n — 1 and 2n 
and 2S — p singlet bonds between sites 2n and 2n-|- 1, for 
every value of n, is a ground state of such a Hamiltonian. 
In terms of the variables u and v, such a state can be 
written as^S 

oc 

*(P) = n [("2n-lW2« - W2«-lU2n)'' 

n— — OC 

("2nW2ri+l " V2nU2n+lT^^^] ■ (30) 

Now, each value of p from to 25 corresponds to a 
ground state of the Hamiltonian; hence there are 25* + 1 
ground states. The case 5 = 1/2 corresponds to the MG 
modelii, while the case 5 = 1 corresponds to the model 
studied in Sees. II and HI. The states in Eq. (|5D|) have 
appeared in the literature as variational ground states of 
a dimerized spin-5 chain, with the integer p changing as 
the dimerization parameter is varied^S. 

One can now consider excitations which are 'domain 
walls' interpolating between ground states ^(pi) on the 
left and ^'(^2) on the right, where, for instance, pi > p2. 
A state of this kind is 

*2n(Pl,P2) 

n-l 

= W [(^t2m-l^'2m - W2m-lW2m)^' 

m— — 00 

X (M2mW2m+l " W2mW2m+l)^'^"^'] 

00 

m—n 

X {u2mV2m+l ~ V2mU2m+lf'^^^^] ■ (31) 

This state has S^^^ = {pi ~ p2)/2 due to the factor of 
yPi-P2 a,t site 2n ~ 1. We can now superpose states 
like this to form a momentum eigenstate, and calculate 
its variational energy. A similar procedure can be used 
to construct excited states interpolating between ground 



One can construct spin models in higher than one di- 
mension in which the excited states exhibit spin fraction- 
alization. Two examples are as follows, 
(i) Consider a spin-1 model on a square lattice in which 
the Hamiltonian 7? is a sum over Hamiltonians H\2 of 
squares for which the ground state has at least two sin- 
glet bonds in each square^l; TJp must be a sum of the 
projection operators P3 and P4 for the total spin of a 
square. The ground states of H consist of a number of 
unbroken lines of singlet bonds such that each square 
has exactly two such lines running along two of its sides. 
Each line of singlet bonds can either extend all across 
the system or form a closed loop. In the limit of large 
system size N ^ the number of ground states grows as the 
exponential of \/N . Hence the entropy per site vanishes 
at zero temperature, even though the number of ground 
states goes to infinity in the thermodynamic limit. Next, 
we can consider excited states in which one of the lines 
ends at a free spin 1/2 at one site; this leaves one square 
unsaturated. Two such excitations are shown in Fig. 4. 
One can then consider variational states in which the free 
spin 1/2 is allowed to move around the lattice in order 
to reduce its energy. 



FIG. 4: Spin-1/2 excitations in a square lattice. For PBC 
there has to be an even number of such excitations. The 
crosses mark the unsaturated squares which have less than 
two singlet bonds running along their sides. 



(ii) Next we consider a spin-1 model on a triangular lat- 
tice in which the Hamiltonian is a sum over Hamil- 
tonians H/\ of triangles for which the ground state has 
at least two singlet bonds in each triangle; -/Ja must be 
the projection operator P3 for the total spin of a trian- 
gle. The ground states of H consist of unbroken lines of 
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singlet bonds such that each triangle has exactly one line 
running along one of its sides. Once again, the number 
of ground states grows as the exponential of \/N for a 
system with N sites. There are excited states in which 
one of the lines ends at a free spin 1/2 at one site; this 
leaves one triangle unsaturated. The free spin 1/2 can 
again move around so as to reduce its energy. 

V. CONCLUSIONS 

We have introduced a Hamiltonian for a spin-1 chain 
which has three degenerate ground states, two of the 
MG type and one of the AKLT type. The lowest en- 
ergy excitation carries spin-1/2 and interpolates between 
the AKLT state and one of the MG states; it has a gap 
Aii^ ~ 2.38J. In the thermodynamic limit N ^ oo and 
temperatures much lower than AE/kB, the system will 
consist of a dilute gas of the spin-1/2 excitation a^^i^^ . 
Hence a quantity like the magnetic susceptibility will go 
as X ~ exp{—/3AE) at low temperatures. The spin-1/2 
nature of these excitations can, in principle, be observed 



in ESR experiments. 

Although the model has three ground states, they will 
not appear with equal weights in the limit of very low 
but non-zero temperature. Since the spin-1/2 excitations 
interpolate between the AKLT state and either one of the 
MG states, we expect that half the chain will be in the 
AKLT state, and a quarter will be in each of the two MG 
states. This implies that the structure factor S{q) at very 
low temperatures will be given by [S^ {q) + S^^^{q)]/2, 
where S^{q) and S^^^{q) are given in Eq. 

Finally, we have indicated how the spin-1 chain with 
spin-1/2 excitations can be generalized to both higher 
spins and higher dimensions. This provides one particu- 
lar way of realizing the idea of spin fractionalization. 
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